Pair-supersolid phase in a bilayer system of dipolar lattice bosons 



o 
o 

<N 



C. Trefzger ^, C. Menotti ^, and M. Lewenstein ^'^ 
^ ICFO - Institut de Ciencies Fotoniques, Mediterranean Technology Park, 08860 Castelldefels (Barcelona), Spain 
^ CNR-INFM BEC and Dipartimento di Fisica, Universitd di Trento, 1-38050 Povo, Italy 
^ ICREA and ICFO - Institut de Ciencies Fotoniques, 
Mediterranean Technology Park, 08860 Castelldefels (Barcelona), Spain 

The competition between tunneling and interactions in bosonic lattice models generates a whole 
variety of different quantum phases. While, in the presence of a single species interacting via on-site 
interaction, the phase diagram presents only superfluid or Mott insulating phases, for long-range 
interactions or multiple species, exotic phases such as supersolid (SS) or pair-superfluid (PSF) 
appear. In this work, we show for the first time that the co-existence of effective multiple species 
and long-range interactions leads to the formation of a novel pair-supersolid (PSS) phase, namely 
a supersolid of composites. We propose a possible implementation with dipolar bosons in a bilayer 
two-dimensional optical lattice. 
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The possibility of engineering lattice models with 
ultra-cold gases in optical lattices is considered one of the 
most promising routes in the search for exotic quantum 
phases which escape clean demonstration in condensed 
matter systems. To the aim of the present work, where 
we show the existence of a pair-supersolid (PSS) phase, 
it is particularly important to briefly introduce the su- 
persolid (SS) and the pair-superfluid (PSF) phases. 

The question whether superfluidity and broken trans- 
lational symmetry can coexist, leading to supersolidity, 
has been intriguing theoretical and experimental physi- 
cists for the last 50 years (see e.g. [ll, i, H, 0, [s, 6]). 
Exact Quantum Monte Carlo simulations have demon- 
strated the possibility of a SS ground state in lattices 
0, 0, The experimental realization of a strong 

dipolar Chromium condensate [ll| and the recent pro- 
gresses towards a deg enerate gas of heteronuclear polar 
molecules [12|, [l3|, IHf put cold gases with long-range in- 
teraction in optical lattices among the best candidates 
for the creation of the SS phase [HI, [IB, [l3 • 

The second important issue concerns particle vs pair 
condensation in Bose gases in presence of attractive in- 
teractions [18]. It has been shown [19, ^201,1211, [13,0,12^, 
that bosonic mixtures with inter-species attraction can 
actually sustain a pair-superfluid phase (PSF) without 
either collapsing or phase separating. The optimal can- 
didates for the realization of the PSF phase are bosonic 
binary mixtures in optical lattices with inters2ecies in- 
teractions tunable via a Feshbach resonance [25[ or, al- 
ternatively, bilayer optical lattices of dipolar particles, 
which create an effectively two-species system if tunnel- 
ing between the two layers is suppressed [2^ . The " inter- 
species" interaction is provided by the long-range inter- 
action, which couples the two separate sub-systems. 

However, dipolar gases offer the further opportunity of 
studying the combined effect of long-range interactions 
and inter-species coupling. In this paper, we show for 
the first time that the presence of both intra-layer repul- 
sion and inter-layer attraction allows for a pair-supersolid 
phase (PSS), defined as a supersolid phase of compos- 
ites. This phase, which joins the characteristics of the 
supersolid phase and the pair-superfluid phase, can be 



obtained in a system of dipolar bosons populating two 
decoupled two-dimensional (2D) optical lattice layers. 



In order to demonstrate the existence of the PSS, we 
study the effective Hamiltonian i^eff in the low-energy 
subspace of pairs, using a mean- field Gutz wilier approach 
and exact diagonalization. Moreover, there exist an accu- 
rate mapping at low densities of the effective Hamiltonian 
onto the single component extended Hubbard model f\] 
the fact that the latter model supports SS implies the 
existence of PSS in our case. 

We consider polarized dipolar particles in two decou- 
pled 2D layers (see Fig. [1]). This geometry can be ob- 
tained by using anisotropic optical lattices or superlat- 
tices, which can exponentially suppress tunneling in one 
direction. The in-plane dipolar interaction is isotropic 
and repulsive. The inter-layer interaction depends on 
the relative position between the two dipoles, but is 
dominated by the nearest neighbor attractive interaction 
W < {) between two atoms at the same lattice site in 
different layers. In the following, we include only nearest 




FIG. 1: Schematic representation of two 2D optical lattice lay- 
ers populated with dipolar bosons polarized perpendicularly 
to the lattice plane. The particles feel repulsive on-site U 
and nearest neighbor [/nn interactions. Inter-layer tunneling 
is completely suppressed, while a nearest neighbor inter-layer 
attractive interaction W is present. 



2 



neighbor in-plane (^7nn) and out-of-plane (W) dipolar 
interactions. The relative strength between /7nn and W 
can be tuned by changing the spacing d± between the two 
layers, relative to the 2D optical lattice spacing d. Due to 
the dependence of the dipole-dipole interaction like the 
inverse cubic power of the distance, the ratio |VK|//7nn 
can be tuned over a wide range. While it can be negligi- 
ble for d± ^ d making the system asymptotically similar 
to a single 2D lattice layer, it can also become relevant 
and give rise to interesting physics, not existing in the 
single layer model [2^, [13, S [Sfl, [3l|, El . 

We start from the generalized extended Bose-Hubbard 
Hamiltonian 



(1) 



(ij) 



where a = a^b indicates the two species (which in the spe- 
cific case considered here are atoms in the lower and up- 
per 2D optical lattice layer respectively), U is the on-site 
energy, /7nn the intra-layer nearest neighbors repulsion, 
W the inter-layer attraction, J the intra-layer tunneling 
parameter and /i the chemical potential. The parameters 
U and J are equal for the upper and lower layers and the 
chemical potentials fi are the same, since equal densities 
in the two layers are assumed. The symbols and (ij) 
indicate nearest neighbors. 

The PSS is characterized by vanishing single particle 
order parameters (a) = {b) =0, and non vanishing pair 
order parameter {ab) ^ 0, coexisting with broken trans- 
lational symmetry, namely a modulation of both density 
and order parameter on a scale larger than the one of 
the lattice spacing, analogously to the supersolid phase. 
The physics leading to the formation of composites re- 
lies on second order tunneling and takes place in the 
low-energy subspace where single-particle hopping is sup- 
pressed. The theoretical description of the PSS phase 
cannot be based on standard mean-field theory, which 
accounts for particle hopping through the replacement in 
the Hamiltonian of the single particle creation and de- 
struction operators by their expectation values, because 
in this approximation second order tunneling completly 
vanishes. 

A successful way to account for second order tunnel- 
ing is to write an effective Hamiltonian in the subspace of 
pairs and include tunneling through second order pertur- 
bation theory [1^, [III, [Hi . The validity of the effec- 
tive Hamiltonian relies on the existence of a low-energy 
subspace well separated in energy from the subspace of 
virtual excitations, to which it is coupled via single par- 
ticle hopping. Such second order couplings are related 
to the super-exchange interaction, recently measured in 
a series of double- well systems [34]. 

The low-energy subspace of pairs is spanned by all clas- 
sical distributions of atoms in the lattice |a) = Hi l^^^ ^i) 



with equal occupation of the two species a and b. For 
(/7 + W^jU 0, asymptotically all classical states \ol) 
become stable with respect to single particle-hole excita- 
tions and develop an insulating lobe at finite J. The 
energy of single particle-hole excitations is of the or- 
der of /7 at J = and is given by the width of the 
lobes at finite J (see e.g. thin blue lobes in Fig|2]). 
This situation has to be compared to the single layer 
situation {W = 0) with nearest neighbor interactions, 
where only uniform Mott phases and checkerboard in- 
sulating phases are stable. The relevant virtual sub- 
space is obtained from the states |a) by breaking one 

composite, namely 17^^"^) = (i};^aj\a)I^Jn^(ni + 1) and 

l7?> = ?'l&.l«)/Vn,K + l). 

In the pair-state basis, the matrix elements of the 
Hamiltonian in second order perturbation theory are 
given by 



(alFeffl/J) = {a|i^o|/3> 
-i^{a|Fi|7){7l^fi| 



(2) 



where i^o, given by the interaction terms of the Hamil- 
tonian ([T]), is diagonal on the states and the single 

particle tunneling term Hi = — ^l^j + ^l^j 

treated at second order. 
For a given state 

E^^^ -Ea = U^{U^W){mi- rrij) + U^^Am^^, (3) 

with Am'^^ T^ik), - E(/c), - 1, where 
indicates the pair occupation number at site i. For 
U+W, /7nn <^ the denominators E^.. — E^ are all 
of order [/, which leads to 
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(0) 
eff 



Ho 



2J2 
U 



E 

(ij) 



^mi{mj + 1) + cjcj 



(4) 



^iff by expanding (|3]) at higher 



where q and are the pair destruction and creation 
operators such that cj I mi) = {mi-\-l)\mi) . One can easily 

obtain corrections to 
orders in {U + W)/U and U^n/U 

In this work, we provide a mean-field solution to effec- 
tive Hamiltonian (|4]). We perform a perturbative treat- 
ment at first order in the pair order parameter = (c) 
to determine the boundaries of the insulating lobes. Fur- 
thermore, we solve the time dependent Gutz wilier equa- 
tions in imaginary time to determine the nature of the 
superfluid phases outside the insulating lobes. We point 
out that in spite of its simplicity, the mean-field treat- 
ment of the effective Hamiltonian is able to include im- 
portant features like the re-entrant behavior of the lobes 
(as predicted by exact t-DMRG calculations for the ID 
geometry in [26!\). 
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FIG. 2: Pair insulating lobes for v = 0,1/2,1,3/2 (thick 
lines); Lobes with respect to single particle- hole excitations 
(thin blue lines) for the dominant configurations in the ground 
state at J = 0.05[/ and /x — — 0.4375?7, namely rrii — and 
rrij = 1,2 (for i, j nearest neighboring sites). The inset shows 
a zoom of the pair phase diagram. 



For any non-zero < 0, the lowest possible exci- 
tations on top of a classical configuration of pairs are 
obtained by adding (removing) one pair at site i. The 
energy costs, respectively given by £^p(J) = — 2/i + 
2/7m, + (2m, + l)W + 2V^{^ - {2J^ /U) E{k)S^^k + 1) 
and Ei^{J) = 2/i - 2U{mi - 1) - {2mi - 1)W - 2V^^^ + 
(2J^/[/) (2m/c + 1), depend quadratically on the 

tunneling coefficient J. In the previous expressions, 
^dip > is the dipole-dipole interaction of one atom 
placed at site i with the rest of the particles belonging 
to the same layer. The mean- field order parameters ipi 
satisfy 



(rrii + 1)^ 



(5) 



with i/ji = X](/c)i being the sum of the nearest- 
neighbour order parameters. Using Eq. ([5j), one can cal- 
culate the mean-field lobes for any given configuration 
of pairs in the lattice. The lobes for the checkerboard 
and doubly occupied checkerboard are shown in Fig. [3] 
for the 0th (full lines) and 1st order (dashed lines) ef- 
fective Hamiltonians. The comparison between the two 
shows that, for the parameters considered here, the 0th 
order already captures the physics accurately. The 
dependence of the energy of the elementary excitations 
is at the origin of the re-entrant behavior of the lobes. 



Based on the Gutzwiller Ansatz for the pair wavefunc- 
tion \^) = Yli Xlrn 1^' predict the existence of 
three different phases: insulating phases, PSF and PSS. 
For U -\-W < zU^^^ being z the number of nearest neigh- 
bors, the insulating phases show checkerboard ordering 
not only at filling factor u = 1/2, but also at filling fac- 
tor V = 1. Outside the insulating lobes, depending on 



density and tunneling, we find either PSF or PSS. The 
stability analysis of the PSS phase against phase sepa- 
ration is beyond our mean- field treatment. However, for 
small U < 2;/7nn, by doping the checkerboard above 
half filling, the extra pair goes on an already occupied 
site and the analogy to the single-species extended Bose- 
Hubbard model [7] suggests that the system stabilizes to 
a PSS phase. Instead, the PSS phase at density lower 
than 1/2 (which our mean- field treatment predicts only 
in a very small region close to the tip of the v = 1/2 
lobe) should be unstable towards phase separation. Nev- 
ertheless, the inclusion of more neighbors to the in-plane 
dipolar interaction is expected to remove this instability 

To get reliable results, one should combine the 
Gutzwiller predictions with an estimate of the limits of 
validity of ^ beyond which the subspace of pairs 
looses its meaning. For each point of the phase dia- 
gram, from the ground state Gutzwiller wavefunction, 
we select the dominant classical configurations with the 



criteria > 0.05 and 111^/^1 > 0.02 |36| and cal- 

culate for each of these configurations, the lobe with re- 
spect to single particle-hole excitations. If the system 
at this given point of the phase diagram turns out to 
be stable against all dominant single particle-hole ex- 
citations (in other words, if this point is inside all se- 
lected single particle- hole lobes), is considered valid 
at that point. This procedure is shown for J = 0.05/7 and 
11 = -0.4375/7 in Fig. [1 In Fig. H we show the result- 
ing phase diagram for /7nn = 0.025/7 and W = —0.95/7. 
The shaded area represents the region of PSS, compatible 
with the above validity conditions for H^^ . Interestingly, 
the ^-like shape of the allowed pairing regions (PSF and 
PSS) matches the ones found in ^26i] . 



The existence of the PSS phase is supported also by 
the exact diagonalization of H^^ for up to 6 pairs on 
a 2 X 8 lattice with periodic boundary conditions. At 
small J, the ground state is almost doubly degenerate. 
This can correspond to two insulating checkerboards (for 
half filling) or two checkerboard supersolids shifted by 
one lattice constant, with a large gap to excited states. 
This double degeneracy, together with a finite measure 
of the coherence, provided by a non vanishing expecta- 
tion value of the tunneling term, are indications of the 
existence of the PSS. At larger J, the quasi-degeneracy 
of the two lowest eigenstates disappears and the ground 
state becomes well separated in energy from all excited 
states, indicating the crossover to the PSF. It is impor- 
tant to remark that in the absence of nearest-neighbor 
interaction, the ground state is always non degenerate 
and no signature for the PSS is ever found. 

Finally, we observe that, upon appropriate renaming of 
the parameters, H^^ (|4|) can be mapped onto the Hamil- 
tonian used in to demonstrate the existence of the 
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FIG. 3: Phase diagram of the effective Hamiltonian, with 
= 0.025[/, W = —0.95U, which can be obtained for 
d± — 0.37d. The black full lines are the semi- analytic solution 
of Eq.(|5]) indicating the boundaries of the insulating lobes for 
the checkerboard (z/ = 1/2) and the doubly occupied checker- 
board (ly = 1). The black dashed lines are the boundaries 
of the insulating lobes for 1st order expansion of i^eff- The 
shaded area is the PSS phase predicted by the Gutzwiller ap- 
proach. The red line indicates the estimated limit of validity 
of . The blue dash-dotted line indicates the upper limit 
of the region where the mapping onto [7] is almost exact. 

SS phase for soft core bosons with nearest neighbor in- 
teractions. Due to the different action of the tunneling 
terms, namely (n^ + 1, rij — l|ajaj|ni, rij) = (rii + l)nj 
in and (m^ + l,mj — l|cjcj Im^, mj) = (m^ + in 



our case, the mapping of the two Hamiltonians is exact 
when only the number states m = 0, 1 are populated. 
Under this condition, the results of [7] translate to the 
existence of the PSS in our problem. The region below 
the dash-dotted blue line in Fig. [3l which corresponds to 
l/d*^P + > 0.9, Vz, namely a close to exact map- 

ping, includes part of the PSS region. 

Summarizing, we have studied the phase diagram of a 
bilayer system of 2D dipolar lattice gases, in the limit of 
close layers, and demonstrated the existence of a novel 
PSS phase, namely a super solid phase of pairs. The ex- 
istence of the PSS phase has been previously discussed 
for anisotropic t-J models [ST^, but no evidence of it has 
been found. However, the Hamiltonian we discuss in the 
present work differs from the anisotropic i-J spin Hamil- 
tonian in three crucial respects, all of which should favor 
the existence of the PSS phase: (i) it deals with soft- 
core bosons (vs. hard-core); (ii) it considers on-site inter- 
species attraction (vs. nearest-neighbor inter-species at- 
traction); (iii) it includes nearest-neighbor intra-species 
repulsion. For these reasons, we believe that the exis- 
tence of the PSS phase will be confirmed by exact calcu- 
lations, also beyond the limits of validity of our effective 
mean-field approach. 
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